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Abstract 



We solve numerically various types of the gap equations developed in 
the relativistic BCS and generalized BCS framework, presented in part I of 
this paper. We apply the method for not only the usual solid metal but also 
other physical systems by using homogeneous fermion gas approximation. We 
examine the relativistic effects on the thermal properties and the Meissner 
effect of the BCS and generalized BCS superconductivity of various cases. 

I. INTRODUCTION 

It is beyond doubt that the relativistic effects are important in condensed matter physics: 
They are present in various systems. But until recently, not so much investigation for 
relativistic effects in the phenomena of superconductivity were performed, except for taking 
into account the spin-orbit coupling [1]. Recently, the papers of Capelle et al. appeared, in 
which the necessity of the relativistic treatment for superconductivity in condensed matter 
was asserted [2~4]. In the study of the superconductivity under strong relativistic effects, 
they argued the necessity of the full relativistic treatment, beyond that of making small 
corrections like the spin-orbit coupling. They introduced the Dirac-type Bogoliubov-de 
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Gennes theory and discussed some new aspects of the relativistic superconductivity. Based 
upon their works, we try to generahze their methods in this paper. 

In part I of this paper, we developed the formahsm of the BCS and generahzed BCS 
superconductivity in relativistic quantum field theory [5]. By using the Gor'kov formalism, 
we derived various types of gap equations. We also gave the group-theoretical considera- 
tion of the pairing gap functions. We treated not only the spin singlet pairing states, but 
also the spin triplet finite angular momentum pairing states by introducing the relativistic 
generahzed-BCS scheme. 

We try now to complete our program and work out the numerical study of the gap 
equations for various cases, by using the formalism developed in part I. To examine the 
relativistic effects on the gap function, thermodynamics and electromagnetic response, we 
have to study various cases by performing numerical calculations. In Sec. II, we examine 
the excitation energy spectra of quasiparticles in the relativistic theory. This examination 
will be useful for investigation and interpretation of the relativistic gap equations and the 
thermodynamics of the superconductor. In Sec. Ill, we solve various types of gap equations, 
given in part I of this paper. We consider not only the situation of the usual solid metal, 
but also other situations where the system has a larger Fermi energy. By studying these 
systems, we depict the strength of the relativistic effects on the gap equations under an unified 
perspective. Because we choose the theoretical framework as the BCS and generalized BCS 
theory, we can materialize it. In Sec. IV, the thermodynamics and the Meissner effect in 
the relativistic BCS and generalized BCS superconductivity are discussed. Finally, in Sec. 
V, we give a summary of the present investigation and provide further possible extensions 
to be studied for the pairing properties of matter. 
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II. THE EXCITATION ENERGY SPECTRA OF QUASIPARTICLES AND THE 

DENSITY OF STATES 



In this section, we examine the excitation energy spectra of the quasiparticles, and the 
density of states (DOS). Especially we concentrate on the scalar and vector pairings, A*^ and 
Aq [5]. As discussed in part I [5], the pairing gap function A4X4 = (V'V'^) is decomposed 
by 16-dimensional complete set of 7 matrices. After constructing the gap equation for each 
component of them, we found the fact that, only A*^, A^, A"^ ( the spacelike components of 
axial vector pairing ) and ^Ja) ( axial- vector-like components of the two-rank antisym- 
metric tensor pairing ) can have nontrivial solutions. Under the treatment given in part I, in 
all the cases of A"^ and A^^^ , the energy spectra of the quasiparticles of these states take the 
similar forms with that of A'^ or Aj', except the differences about the angular dependences 
of the gap functions in momentum space. Thus, we regard the energy spectra of A*^ and A^ 
as having the typical character of the quasiparticles in the relativistic BCS and generalized 
BCS theory. 

First, we give the excitation energy spectra of the quasiparticles of the A'^ and A^ states 
schematically in Fig. 1. £'+(£'_) is the branch of the quasiparticles coming from positive ( 
negative ) energy states. The explicit forms are given as 

E± = ^(Vfc' + m2^//)2+|A5|2, (1) 

for the scalar pairing [2,5], and 

E± = e + + 11^ + \A^\^ T 2^{e + w?)ii' + |A^|2fe^ (2) 

for the 0th component of the vector pairing [5]. We consider first the spectrum In 
the case of A'^, the minimum is found at kmin — kp = V f^'^ ~ ( Throughout this paper, 
we take the approximation to set 11 = ep. This is the case for zero temperature. We 
completely neglect the temperature dependence of // and use it as a simple parameter to 
give the particle density of a system. ), and the energy gap is given by the width 2|A'^|, 
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same as the usual nonrelativistic case. On the other hand, in the case of A^, the minimum 
is found at kmin — \/ f^^ + l^o'P ~ m?p? /{p? + |Ao'p)(> A.^), and the gap width becomes 
2m\/^\/ ^JJl?^^\^^{< 2|A^|). The gap exists also eX kp-, but particles feel smaller gap at 
kmin- Usually > |A^|, and therefore the gap width becomes ~ 2^|A^|. For < 1 ( 
like the case of usual sohd metal ), the gap width becomes nearly 2|A^|. In this case, the 
spectra of A'^ and A^ almost coincide with each other. It is also the case that, when ii — m 
becomes large, the difference of the spectra between A'^ and A^ becomes large. Puther, it 
is clear from the above discussion, if we treat the massless case, the gap will vanish in the 
A^ state. 

The branch of E_ is located above E+ in both cases. When <^ 1, the relative 
difference between and E_ is large compared with the Fermi energy, ii — m{ which gives 
the characteristic energy scale of a system ). Hence, the contribution of the quasiparticles 
of becomes small or can be neglected. On the other hand, when // — m ~ m ( like the 
relativistic plasma ), the contribution of E_ would become large, and we should take E^ 
into consideration for the pairing properties. 

The density of states ( DOS ) is estimated by 

ps{uj) = Pni^p) j di5{uj - E+) 

= Pn{eF) JdE^^6{u;-E+), (3) 

where Ps{uj) is the DOS of the superconducting state, while Pni^p) is the DOS of the normal 
state at the Fermi energy; Pni^r) — p^/p^'^-rr?, and ^ = y/k'^ + m? — p. Here we only 
take into account the contribution of the branch E^. Fig. 2 schematically depicts the DOS 
for the A*^ and A^ states. In the case of A'^, Ps(cj) approaches to infinity aX u ^ ±|A'^|, 
while in A^ state, ps{oj) becomes infinity at a; ^ ±m\/S^\/ ^ p^ + IAq^^. Therefore, when 
lA"^! = |A^|, the quasiparticles in the A^ state is more easily excited than that of the A'^ 
state. 
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III. NUMERICAL SOLUTIONS OF THE GAP EQUATIONS 

In this section, we discuss the numerical solutions of the gap equations. In the numerical 
calculations, we should not use the approximation which is usually used in the condensed 
matter theory: 



Rather, we have to treat the details of the energy spectra in the integration of the gap 
equation in order to compare the relativistic and nonrelativistic cases. Furthermore, we 
have to introduce a cutoff to make the integral of the gap equation finite due to the use of 
the 5-function expression for the interaction for simplicity [5] . We take the energy cutoff A 
in the following way: 



As discussed in part I, the aim of this work is to study the relativistic BCS and generahzed 
BCS theory. For this aim, we set aside the question of the mechanism of the origin of the 
superconductivity. This is the attitude of the BCS and generahzed BCS theory. Due to this 
attitude, we can study the relativistic superconductivity of various systems in an unified 
approach, from the usual solid to the stellar matter. Therefore, we treat the cutoff in a 
general way. In the electron-phonon interaction, we only treat a thin shell near the Fermi 
surface. Since we skip the question about the origin of the attractive interaction, we can 
extend the width of integration as Tq < A < Tp, based on the discussion of the textbook 
of VoUhardt and Wolfle [6], chapter 3. Here ksTq = l/r ( r; the quasiparticle lifetime of 
the Fermi liquid ), while Tp means the Fermi temperature. We integrate the states of the 
inside and outside of the Fermi sphere symmetrically, with moving the cutoff A in this range. 
Then, the gap equation generally has the cutoff A, coupling constant gi, fermion mass m, 
chemical potential and temperature T as the parameters of the model. 

Our gap equations have five parameters. Therefore, there is a large arbitrariness to choose 
these values. In this work, we mainly investigate the relativistic effects in superconductivities 




(4) 





5 



of several physical systems. To determine the conditions of solving the gap equations, we 
present the typical order of physical quantities of several systems, in table I. In this table, we 
consider 5 examples: (i) The usual solid metal ( metallic elements, alloys, compounds, oxides 
and organic metals ) [7,8], (ii) the metaUic hydrogen [8,9,10], (iii) to depict the strength of 
relativistic effects in heavy elements, we introduce an imaginary system of electron gas, in 
which ep — m is almost same as the kinetic energy of Is-electron of uranium atom, (iv) the 
relativistic plasma [11] and (v) the neutron star [12]. To estimate the Fermi energy of the 
systems, we use the nonrelativistic energy-momentum relation ep — m + ior the cases of 
(i), (ii) and (v), while we use the relativistic energy-momentum relation ep = ^kp + for 
the cases of (iii) and (iv). For (ii), the value of = ^ ( Tq; mean particle distance, qb] the 
Bohr radius ) takes the range of 0.5 < < 1.5 [9,10]. Based on the paper of Ichimaru [10], 
the electron system is relativistic at < 0.1 with kBT/^ep — m) <C 0.1 ( at = 0.1, 
ep — m — 5 X lO^eV and ~ 10~^ ). For (iii), we take the mean radius of Is-electron 

as {r)is — 715 x lO^^^cm [13], and estimate the particle density by (|7r(r)5^)~^. In this 
case, we obtain = 0.013. It is known that the relativistic effect ( kinetic energy effect ) 
enhances the screening, and the screening length takes a finite value at rg — > 0, while the 
nonrelativistic Thomas- Fermi length diverges in the same limit [10,14]. Therefore, we shall 
consider the superconductivity in such high density matter, with an assumption of existence 
of some effective attractive interaction. 

Consulting with table I, we put the values of ep to n in the gap equations. Commonly, 
the magnitude of the gap is much smaller than /j, — m. For example, in usual solid metal, 
the ratio |A(r — 0)|/(// — m) becomes 10~^ ~ 10~^. Hence, we search the values of our 
parameters, especially gi and A, to satisfy this condition in solving the gap equations. After 
choosing the values of the parameters, the integration in the gap equation is performed, 
and search the self-consistency condition under the variation with respect to the amplitude 
of the gap. For the integration, we used the numerical package Mathematica version 4.1. 
Usually, this package gives high numerical accuracy. 

First, we treat the spin singlet pairing states. We take the cases of A'^ ( scalar ): 
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E± = ^(£=F/i)2 + |ASp, (7) 
( 0th component of vector ) : 



1 = 



r deeV^^^Ul - , ' -"^ }4^tanh 



+{1+ , ^ =}^tanh^£;_), (8) 

/eV + |A^P(e2-m2) 2£;_ 2 ^ 



= ^Je^ + pi^ + |A^|2 IF 2^eV + |AS'|2(e2 - m^), (9) 



"no-sea" ( neglect the contribution from the negative energy states in A'^ ): 

/ deeVe^ - m^—— tanh ^E^, (10) 



47r2 



= V^(6-/X)2+|A^P, (11) 

and "nonrelativistic" ( substitute y/^~+m^ + ^ "no-sea" ) : 

1 -|?.r>^V2-^(--)^tanhfE, (12) 



E = ^(e - /^)2 + |A^«|2. (13) 

Therefore, we treat these 4-cases of the gap equations for comparison. 

Fig. 3 gives the gap functions at T = as a function of the couphng constant Qq. Here we 
take the cutoff A such that the integration is performed over the inside and the outside of the 
Fermi sphere symmetrically; A/(/i — m) = 1. In the case of the solid, the solutions obtained 
from these 4-types of gap equations discussed above coincide with each other very well. In 
this case, the relativistic effect is negligible at least in the gap function. There is no contribu- 
tion of the quasiparticles coming from the negative energy states. In fact, the gap equations 
of A-^ and A^, Eqs. (6) and (8) almost coincide in the hmit < 1 with |A| < On the 
other hand, in the case of the relativistic plasma, there are large differences between these 
4 solutions. A"^ gives the largest gap, while A^ is strongly suppressed and becomes the 
smallest gap. The gap of the "nonrelativistic" case is also largely reduced compared with 



A"^. The difference of A"^ and "no-sea" is also clear. At go — 15.8MeV~^, the relative ratios 
are 105 for A'^/ "nonrelativsitic" and 81.4 for "no-sea" / "nonrelativistic" . Therefore, in the 
relativistic plasma, we have to introduce the relativistic treatment, and cannot neglect the 
contribution of the negative energy states. For the case of the neutron star, the solutions 
of the 4-cases are close from each other. But we observe small differences between them. 
At go = 3.55 x 10~^MeV~^, the relative ratios are 1.16 for A'^/ "nonrelativistic" , and 1.15 
for "no-sea" / "nonrelativistic" . This indicates that the contribution of the negative energy 
states to the gap function is less than 1 percent in this case. It is clear from Fig. 3, the 
solutions reflect the nonperturbative effect and behave almost e~^^^° in all cases. 

In Fig. 4, we show the dependence of the cutoff A, namely the integration width, in 
the solutions of 4-cases at T = 0. Here, x = 2A. We squeeze the width to a small region 
containing Fermi surface. At A/(/i — m) = 1, we integrate all the inner region of the Fermi 
sphere and outer region symmetrically. Same as in Fig. 3, the solutions of the 4-cases 
coincide very well in the case of the solid, and |A| <C A is always retained. In the cases of 
the relativistic plasma, when we squeeze the cutoff, the gaps of A*^ and "no-sea" come close 
to each other. But the relative ratio between A'^ and "nonrelativistic" still remains large. 
We also observe that in the case of the neutron star, the solutions of the 4-cases come close 
to each other when we squeeze the cutoff. In all cases, the solutions almost linearly depend 
on X. Therefore, we conclude that the solutions behave almost |A(T' = 0)| oc Ae~^/^° ( like 
the nonrelativistic BCS theory [8,15,16] ) in 4-cases of various systems. 

In Figs. 3 and 4, the values of |A(T = 0)| are given almost in the same range. From the 
results of Figs. 3 and 4, we find the values of go and A to provide reasonable magnitude of 
the gap, and then we solve the T-dependent gap equations. Fig. 5 gives the temperature T 
dependences of the solutions in the 4-cases. For the solid, the solutions of the 4-cases almost 
coincide. The critical temperature Tc of Fig. 5(a) is 20.7 K. The case of metallic hydrogen 
( here we take = 0.5 ), we find a tiny difference between A*^ and "nonrelativistic" cases, 
and Tc ~750 K. In the case of the relativistic plasma, Aq and "nonrelativistic" are largely 
reduced. The difference of A"^ and "no-sea" is also large. Tc of A"^ is 0.0820 MeV = 9.5x10^ 
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K. This Tc is 0.16 times the electron rest mass. For the neutron star, Tc of Fig. 5 for A'^ is 
1.86 MeV, while 1.71 MeV for the nonrelativistic case; thus the magnitude of the relativistic 
effect is less than 10 percent. 

Next, to describe our numerical results of the singlet states of various electron systems 
in an unified manner, we solve the gap equations with very fine variation of the chemical 
potential. We show the results in Fig. 6. The value of qq is determined so as to provide the 
ratio I A(T = 0)|/ {ji — m) — 10~^ in the solution of the "nonrelativistic" case obtained from 
Eq. (12). The cutoff is always set as A/ [ji — m) — 1. By using these values for go , we solve 
the gap equations (6), (8), (10) and (12), and obtain the values of |A(T = 0)| and Tc. The 
deviations from "nonrelativistic" give the relativistic effects. The relativistic effects become 
significant for — m ~ lOkeV in both quantities. At ji — m — O.SMeV, Tc for A'^ case 
becomes 0.077 times the electron rest mass. We obtain an important result that, by using 
the data in Fig. 6, the BCS universal constant |A(T = Q)\/kBTc^l.l& [8,15,16] is almost 
completely filled in A'^, "no-sea" and "nonrelativistic" in all ji — m, while discrepancies 
are obtained aX fi — m > 50keV for A^ state. Thus, we conclude that A*^, "no-sea" and 
"nonrelativistic" cases obey the BCS-like temperature dependence in all situations. 

From the above results, we conclude that if the chemical potential is large enough, 
generally the gaps of the A*^ case and "no-sea" are larger than that of the nonrelativistic 
case, and Tc also becomes higher. Especially in the case when the attractive interaction 
requires the treatment covering the wide range of inner and outer regions of the Fermi 
sphere, those relative differences become significant. It is well known that, the DOS of the 
relativistic theory in normal state eV — m? is larger than that of the nonrelativistic one 



yj2m^{e — m). This fact also contributes to enhance the solutions of the relativistic gap 
equations as compared with nonrelativistic counterpart. The relative difference between the 
"no-sea" and "nonrelativistic" cases has its origin in this effect. We expect from the relative 
relation among the magnitudes of the gap, the condensation energy, the lowering of the free 
energy and the thermodynamic critical magnetic field become larger in the relativistic case 
than in the nonrelativistic case. When (// — m)/m <S 1 ( like the usual solid metal ), we 
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cannot distinguish A'^ and A^. In such a case, the presence of the state is "hidden". 
But actually, A"^ and can be mixed under 0(3) rotational symmetry ( as discussed in 
part I of this paper ). 

The reason of the relation | A'^| > | A^| in all /i is understood in the following way. Based 
on Sec. Ill of part I, each mean field is given as follows: 

A^oc(^^(-C75W, (14) 
Ao^oc(V^^(-C75)7oV')- (15) 

Here, T means the transposition of a matrix. In the standard representation, is given as 



(16) 



where ip is the large component, while % is the small component. Thus, from the relations 
of (14) and (15), we obtain 

A^ = ip^{i(T2)^ + X^{i<y2)x, (17) 

Al = ^'^{ia2)v-X^{i(T2)x- (18) 

Here i(72 expresses the spin singlet pairings clearly. At ^ m, the small component vanishes, 
and A"^ and A^ coincide. On the other hand, when — m ~ m and the kinetic energy 
becomes large, </? and x approach to the same limit ( when m/ ^ — > 0, they coinside with 
each other ), then the difference between the magnitudes of the A'^ and A^ becomes large. 
( We can choose the phase 6^ = in A"^ = \A^\e''^ and A^ = |A^|e'^. ) 

In this work, we do not present the results of the massless case and/or ultrarelativistic 
limit. We mention, however, for the case of Aj^, we obtain an unphysical result at m = 0; 
two solutions appear in our numerical calculation. 

Next we give the results for the spin triplet gap equations in Fig. 7. We solve our 
gap equations, given in (114)~(117) with table I of part I, under the 3 conditions: (a) 
The usual solid metal, (b) the relativistic plasma and (c) the neutron star. In all these 
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cases, at T = 0, the largest solutions are obtained for the ■0ooo state of ^Ja)- This state is 
regarded as the Balian-Werthamer (BW) state [6,8,17] in our theory, as discussed in part I 
of this paper [5]. Thus our theory also predicts that the BW state is the most important 
state, like the nonrelativistic theory. The ■0ioi state of A"^ ( the Anderson- Brinkman- Morel 
(ABM) state [6,8,18] in our theory [5] ) gives smaller solutions than that of the BW state 
at T = 0. The of the BW, ABM and ^(V^ni - V^tn) of coincide very well. In 
all 3 conditions given above ( (a), (b) and (c) ), the values of |A(T = 0)|/A;bT'c are 6.25 
for the BW, while 5.87 for the ABM. Especially for the case of A^^V'ooo and A^V'ooo, their 
relations to each other in the gap equations are the same as that for A"^ and A^ except for 
the coefficients. Therefore, to interprete the results of the spin triplet gap equations, we can 
use the main results of the spin singlet states. We can conclude that, when (/x — m)/m <^ 1, 
the solutions of A"^ and A^^ in the same ipjmjX will coincide very well. On the other hand, 
when fj, — m m, the cases of the transversal states of A^ and the longitudinal states of 
A^^ remain important states; the solutions of other cases given in part I will be strongly 
suppressed and become unimportant. Thus we recognize the fact that, in the relativistic 
theory, the relative orientation between the spin vector and momentum vector ( it is given 
by the helicity of the Cooper pair, as discussed in part I ) determine the magnitude of 
the stabilization of the superconductivity. Under the variation of parameters { gi, A and 
fjL )of the gap equations, the transversal states of A^ and the longitudinal states of A^^ 
behave hke A'^, while other cases given in part I behave hke . It is also the case that, 
in the cases of j = 2 pairings, the solution of their states becomes relatively much smaller 
than those of the BW and ABM states. These results come from the fact that the basis 
i^jmjxio — 2) for expanding these states has components which do not couple with p-wave 
interaction Qi^Yi^Y^^, because of the orthogonality coming from the azimuthal angular 
dependence. If we extend our theory to include f-wave components in the pairing interaction, 
these situations may become different. 

To summarize the numerical results of this section, we would like to mention that it is 
impossible to remove the arbitrariness of the choice of the values of the model parameters 
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completely. Especially, the magnitudes of the solutions of the gap equations are sensitive 
to the coupling constants and the cutoff parameters. Therefore in our work, we can only 
discuss the characteristic features of the physics of the relativistic BCS superconductivity. 
We assert, however, that the various qualitative features we have obtained are true and we 
have to keep in mind the importance of the relativistic effects in some situations discussed 
above. In such situations, the results obtained here should be considered seriously. 

IV. THE THERMODYNAMICS AND THE MEISSNER EFFECT 

In this section, we discuss the thermodynamics, especially the specific heat and spin 
paramagnetic susceptibility, and the Meissner effect. First, we consider the specific heat of 
the superconducting state. It is well-known in the nonrelativistic theory [6,8,19] that the 
temperature dependence of the specific heat C at T <^ is determined by the DOS at the 
Fermi surface. The DOS of the Fermi surface is refiected by the node structure of the gap 
function. If the gap has no node, the specific heat depends exponentially on T ( like the 
usual BCS or the BW state ), while C oc T'^ for the case when the gap vanishes at a fine ( 
the polar state ) and C (x for the case when the gap vanishes at two points ( the ABM 
state ). To estimate the T-dependence of the specific heat in the relativistic theory, we use 
the next formula: 

^ °^ / (0^^ f ^+ + ^- f • ^'^^ 

( for all the details of the derivation, see appendix A ) Thus, C is determined by DOS, T, 
E^, E_ and the derivatives of the Fermi distribution functions. The second term is usually 
neghgible. When [ji — m) / m <^ 1, the excitation spectra of A'^ and almost coincide. On 
the other hand, when ji becomes large, the difference between the excitation spectra of A'^ 
and A^ becomes significant, as discussed in Sec. II. In such a case, the quasipaticles are more 
easily excited in the A^ pairing than in the A*^ pairing. But in this situation, A^ becomes 
unimportant, because in such a case, the solution of the gap equation of the A^ becomes 
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small compared with A'^, "no sea", "nomelativistic" , as discussed in Sec. III. Therefore, the 
effect of the difference in the excitation spectra of the radial direction in fc-space ( like the 
difference of in A*^ and A^ ) does not affect the thermodynamics seriously. We conclude 
that, in the relativistic theory, we can also estimate the T-dependence of the specific heat 
by the node structure of the gap. Therefore, A-^, , BW and ^(Y^ni — ^iln) of A"^ of 
our theory have the exponential dependence, while the ABM state of our theory has T^- 
dependence. In our theory, there is no state which has the polar state like node structure ( 
see, part I of this paper ) . 

Next, we consider the spin paramagnetic susceptibility. Here we only consider A*^, A^ 
and "nonrelativistic" cases; the spin singlet states. We introduce the Zeeman energy for a 
weak homogeneous external magnetic field, and employ the usual treatment, obtaining the 
magnetization 



M 



(20) 



Prom this, the spin paramagnetic susceptibility becomes 



^ o f d^k ,df{E+) df(E_)^ 

Here the second term is negligible. When we neglect this term, the difference of A'^ and 
"nonrelativistic" is only the DOS. Prom this observation, we suppose there is no qualitative 
difference in the T-dependences in A"^ and "nonrelativistic" cases. When — m ~ m, the 
DOS enhances Xs of A*^, but the gap also becomes large compared with that of the "non- 
relativistic" pairing, and thermal excitation becomes more difficult in A'^ than in "nonrela- 
tivistic". In this situation, |A^| is largely reduced and the A^ state becomes unimportant 
( as discussed in Sec. HI ). To retain the comparison of A*^, A^ and "nonrelativistic" 
meaningful, we have to take // — m small enough or a moderate value. In such situations ( 
the usual solid metal, metallic hydrogen and neutron star ), the DOS and E^ of them are 
almost same values. In fact, we obtain numerical results that, in such situations, A"^, A^ 
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and "nonrelativistic" cases give the same T-dependences for Xs-, and behave hke the usual 
Yosida function [8,20]. We also confirm in our numerical calculation, at large [x ( relativistic 
plasma ), the Xs of A*^ also behaves like the Yosida function ( at large — m, it is needless 
to consider the state ), and the second term in (21) is completely neghgible. 

Strictly speaking, the treatment given above is not completely a relativistic one, because 
the Zeeman energy is a first-order term in 1/c ( here c is the velocity of light ). Here, we only 
want to see the effect of the magnetic field in the simple method, which is usually used in the 
nonrelativistic theory. To compensate this incompleteness, we calculate the electromagnetic 
response function, and examine the Meissner effect, especially for the A'^ and A^ states. 
The derivation of the response function is given in appendix B. We find, at a small or a 
moderate value of — m ( the usual solid metal, metallic hydrogen and neutron star ), the 
T-dependence of the numerical integrals of (B17) become almost 1 — (T/T^)^ in the A*^ and 
Aj' states, similar to the behavior of the response in the nonrelativistic theory [16]. This 
behavior is found also for the case of large n — m { the relativistic plasma ) for A*^ pairing. 

In summary, we conclude that, in the T-dependence of several thermodynamic quantities 
and Meissner effect, there is no qualitative difference between the relativistic ECS theory 
and nonrelativistic BCS theory in homogeneous systems. Therefore, we can safely use the 
basic concepts they are used in conventional nonrelativistic BCS and generalized BCS theory 
of superconductivity, for qualitative understanding of several thermodynamics and Meissner 
effect in various systems, from the usual solid metal to the stellar matter. This fact is not 
clear before our work is done. 

V. CONCLUDING REMARKS 

In part I and part II of this paper, we have performed the investigation of the BCS and 
generalized BCS superconductivity in relativistic quantum field theory. After the prepa- 
ration of the Gor'kov equation and the group theoretical consideration of the pairing gap 
functions, we have solved the Gor'kov equations under the conditions which are assumed 
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to have specific types of the pairing gap functions. After that, we have constructed various 
types of gap equations, and have solved them numerically We have discussed the results in 
detail. The thermodynamic functions, especially the specific heat and the spin susceptibility 
were discussed. Finally we have examined the Meissner effect in our theory. We have ob- 
tained and examined various characteristic features of the relativistic BCS and generalized 
superconductivity throughout this paper. 

In part II, we have mainly investigated the kinetic energy effect of the homogenous 
systems in the relativistic superconductivity. This effect clearly appears in the gap equations 
of several situations. On the other hand, the works of Capelle et al. [2~4] treated the systems 
under deep periodic potentials. They considered the inhomogeneous effect in the relativistic 
superconductivity, and new physics they have obtained are based essentially on this effect. 

Now we discuss further possible extensions. In the case of the vector or tensor ^Ja) 
pairings, we have treated only the unitary cases. These treatments should be extended 
to study the case of nonunitary states, to solve the Gor'kov equations and to obtain gap 
equations. In the case of nonunitary state, time reversal symmetry is broken spontaneously. 
We are interested in the several physical properties of the nonunitary states, as recently 
discussed in a nonrelativistic theory [21]. 

The investigation of the collective modes is also interesting in our theory. We have 
introduced new mean fields for pairings, and therefore they may have collective modes which 
have not been known yet. For example, various collective modes in the BW and ABM states 
were studied in detail in nonrelativistic theory [6,8,22~24]. It is interesting to study how 
these modes are modified under the relativistic effects. 

Extension of our theory to the Ehashberg formalism ( strong coupling theory ) [25] is an 
important theme. We should study how the retardation of the pairing interaction modifies 
our results. In the context of the relativistic theory, inclusion of the Coulomb interaction ( 
and/or the photon- mediated interaction ) is also interesting. 

We have obtained the Green's functions of the relativistic superconductivity in this paper, 
with which we gain the abilities to investigate various physical properties. The problems 
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given above may be treated by our Green's functions. These investigations demand us of 
preparation of the next stage, which is outside of this paper. 
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APPENDIX A: DERIVATION OF THE SPECIFIC HEAT 

We write down the thermodynamic potential ( in unit volume ) in our theory: 

Q,(r) = Q„(T = 0) + £;c(T = 0) 

-^BTtr / ^ln(l + e-^^+)(l + e-'^^-), (Al) 

where the first term gives the normal state thermodynamic potential at T = 0, the second 
term gives the condensation energy at T = 0, and the third term gives the contribution of 
the Bogoliubov quasiparticles. The entropy is given by 

/d k c 
^{(1 - /(£;+)) ln(l - f{E^)) + f{E^) ]nf{E^) 

+ (1 - f{E.)) ln(l - f{E_)) + f{E_) ln/(E_)}. (A2) 

This is the formula for the entropy of quasiparticle ideal gas. Here, f{E±) is the Fermi 
distribution function for the Bogoliubov quasiparticles: 

^(^*) = exp(A) + r (^'' 

At A — > 0, Eq. (A2) gives the entropy for the normal state. This indicates a second-order 
phase transition. We yield the formula for the specific heat in the relativistic theory: 
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dT 



d/3 



tr 



dE 



f d'k , r,^ odE+,^ ,2/3^ 



dE. 



+ iE-+P^-^)E-Sech^^E.}. 



(A4) 



Futhermore, at T <S Tc, the gap function satisfies the condition A(T) fa A(T = 0), and 
therefore dE±/dj3 = 0. We intend to study the temperature dependence of the heat capacity 
at T ^ Tc- Hence, we use a formula which neglects the temperature dependence of the gap 
function. Therefore, we obtain (19) in Sec. IV. 

APPENDIX B: RESPONSE FUNCTION FOR THE MEISSNER EFFECT 



Our starting point is the Gor'kov equation under the presence of an external field, (13) 
of part I [16,26,27]. In this equation, we regard as a perturbation, and see the variation 
in the one-particle propagator in the first order. We decompose Green's functions as Sp — 
Sf^ + S^p\ F = F(0) + F(i), F = + and put them into the Gor'kov equation, 
retaining only the first order terms. Here (0) indicates the 0th order, the part which satisfies 
the Gor'kov equation with no external field, while (1) indicates the variation with respect 
to the external field. Due to the choice of the Coulomb gauge, we assume that the mean 
field is rigid under the external field [26] . We obtain 



v 



i'l^dn — m + ^^ji A 

A i'y^^d^ + m — 7°^/x 



S^^\x,y) -iF(^\x,y) 



\ 



V 



/ 



sP{x,y) -iF^'\x,y) 



-eYAi,(x) 
e-f'^^Ai.ix 

and write it in the form of the Dyson-type integral equation. 



-iF(^\x,y) -S^F\y,x 



(Bl) 
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X 



^ -iF^^\x,z) 



-iF^'^\x,z 



Sp (z, x) ^ 



-eYA^{z) 
e^^'^A^iz) 



I 



St\z,y) -iF^^\z,y) 



\ 



-tF('Hz,y) -SPiy^zf 

(B2) 



We take the first order variation, for example for Sp- 



sP{x,y) = -e/ d'z{sP{x,z)rA^{z)SP{z,y) + F('\x,z)r^A,{z)F('\z,y)). 



(B3) 



The induced current now becomes 



j^ix) ^ie hm tr-fi^sP {x,y) 

= -ie" J d^ztr^^^'SPix, z)YSP{z, x) + ^^F^''\x, z)-f^^ F^^'^z, x))A,{z) 
^e^ j d^z\l^,{x,z)A%z). (B4) 



We obtain the polarization function as 



Tl^.{x,z) = -ieHv(ji'SP{x,z)YSP{z,x)+^i'F^''\x,z)Y'^m{z,x)). (B5) 



This is given in the Fourier transform as 



= -iehr J ^,{rsP{k)YSP{k + q)+rF^'\k)Y^F('\k + q)). 



(B6) 



Following the same way, we obtain formulae in the Matsubara formalism. The induced 
current becomes 



j'^{x) = -e lim tr^'^S^^\x,y). 



(B7) 



The polarization becomes 



(B8) 



and in the Fourier transform 

1 r d^k 



(B9) 



Hereafter we use only the Matsubara formahsm. 



In the gauge invariant theory, the transversal condition for the polarization, q^U^,^ = 0, 
comes from the current conservation. In the case of zero temperature quantum field theory 
for the vacuum, from the Lorentz invariance, it takes the form n^,^(q') = {Qfiuff — Q.tiQ.v)^{(f) i 
which represents the 4-dimensinal transversality In our case, the system is under finite- 
temperature and finite-density. We take a specific coordinate for the imaginary time and 
select a suitable rest frame for the description of the system. Then in our case, the Lorentz 
invariance is lost. Only 0(3) rotational invariance in fc-space remains. Then we use the 
following projection operator ( here q= {qo = i^^^TT = ia;„, q) ) [28]: 

= ^i^H5)n(^Hgo, \q\) + PiP(g)n(^)(go, \q\l (BIO) 

and from P(P + P(^) 

— Qixv ~ Q.ijl(Iv / (f' 1 we derive a 3-dimensional transversal condition, and 



then 



q^q^ 



( 



\ 



(Bll) 



(B12) 



^ \ Qijq^ + qiqj j 

P^P is for the 3-dimensional transversal condition, while P^^^ is the 3-dimensional longitu- 
dinal condition, and each of them satisfies q^H^^j^") — 0. Using them, we obtain 



n(^n?o,kl) = -Voo(?), 
q 

n(^Hgo,kl) = knj:(g) + ^n, 



(B13) 
(B14) 
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Next we choose the Coulomb gauge q^Ai{q) — and also take Aq — 0, and the response 
becomes only 3-dimensinal transversal: 

Mq) = -n^^\q)A{q). (B15) 

For the Meissner effect, we only consider the static field qo — O, then 

j(0, q) = -^(Hn + + n33)(0, q)A(0, q) 

= -n(0,q)A(0,q). (B16) 

We calculate this response by our Green's functions. We concentrate on the case of the 
scalar A'^ and the 0th component of the vector A^. For the finite q, it becomes rather 
complicated calculation. Here we only want to study the Meissner effect, then we take the 
London limit q ^ 0. After the frequency summation, we obtain 

d^p 1 



/dp L 
j2^nEl{p)-El{p)Y 



^ ;N{E+{p),p)^sech'^E+ip) 



E4py ' ^2 2 

h' 

1 5E4pf-E4p)\^^^ „^.. 



+ iP) , P) f sech^ ^E_ ip) 



where N{po,p) for the scalar A'^ case is 

+ (Po + + "^^)(5p^ + 9m^) 

-(Po + 6a*Vo + A*^ + (P' + m^) (P' + 3m2) 

+ \^'\'[-{pI - fi'){3pl + 9//^ - 2p' - 6m') + {p' + m'){p' - 3m')] 

+ |A^|4[-3(p2 - 3//) + 5p' + 3m'] + 3|A^|6}, (B18) 



and for the vector A^ case is 
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+ (Po + + m2)(5p2 + 9^2) 

-|Ao^n3(3p^ - 2^iy, - I,') - {2>pI - ^i'){p' + 3m') 

+(Po + A*')(P' - 9m') - 3(p' + m')(p' - 3m') + Ap'{p' - m')] 

+ |A^|^[3(3p2 - ij,^) + 5p' - 9m'] - 3| A^|^}. (B19) 

In (B17), the first and second terms become zero at T = 0, and only contribute the neighbor 
of the Fermi surface to the integraL The integral of these two terms always converge. These 
terms have their origin on the terms which describe coUisionless elementary excitations [29] . 
The third and fourth terms give diverged integrals, and we interpret these terms as coming 
from the terms which describe the creation- annihilation of the Cooper pair [29] . The third 
and fourth terms have small temperature dependences. In the finite-temperature field theory, 
is devidcd into U^^ = H™"**^'' + n;;^^""'" with the definition n;;^^""'" = lim^^o,r^o n^^, 
and the usual renormafization procedure is adopted to the nj^"*^"""* [28,30]. Here we also use 
this method to handle the divergence of (B17). Therefore, 

n^.(g) = {U,,iq) - n— (g)} + R— (g), (B20) 

and 

W J {27r)^^'^^k^-m^^''{k + q)^-m^- ^ ' 

It is well-known that, after the method of usual gauge- invariant regularization and renor- 
mafization is performed to nj^"™""*(g), this term becomes {QhvQ' — QhQv)^{q'), and vanishes 
at go = 0, qr ^ 0. Therefore, for our aim, fimg_o{n^^(0, q) - U^^'^O, q)} should be cal- 
culated. By using this function, we examine (B16) at g ^ 0. This examination corresponds 
to treat the London equation at q ^ 0. Thus, the temperature dependence of the integral 
(B17) corresponds to the temperature dependence of the supercurrent in the relativistic 
superconductivity. 
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FIGURES 

FIG. 1. A schematic figure of the excitation energy spectra of the quasiparticles of the A'^ and 
Aq states. and correspond to the quasiparticles coming from the positive energy states 
and the negative energy states, respectively. In E^, the minimum energy appears at a momentum 
denoted by kmin, which is different between the A'^ case and Aq^ case, which are explicitly written 
below Eq. (2) in the text. 

FIG. 2. A schematic figure of the density of states ( DOS ) for the spin singlet pairing states, 
A'^ and Aq , normalized by p„(eir) as a function of energy. The region around the gap is depicted. 
We only take into account the contributions of the quasiparticles coming from positive energy 
states. In the case of A*^, the DOS approaches to infinity at a; — > ztlA*^!, while in Aq^, the DOS 
becomes infinity at a; — ibmlAo'K/x^ + |Ao'P)'~2. 

FIG. 3. go dependence of the pairing gap of spin singlet states at T = 0. (a) The usual solid 
metal. We set jj, — m = 5eV (r^ = 3.2). (b) The relativistic plasma. We set /x — m = 0.5MeV 
{rg = 0.0082). (c) The neutron star. We set fi — m = 20MeV. The masses are 0.5 MeV for electron 
and 1 GeV for neutron. All calculations are performed with the cutoff A/(/i — m) = 1. 

FIG. 4. The cutoff dependence of the pairing gap at T = 0. (a) The usual solid metal. We set 
/X — m = 5eV and go = 9.48 x 10~^eV~^. (b) The relativistic plasma. We set iJ. — m = 0.5MeV and 
go = 23.6MeV-2. (c) The neutron star. We set /x - m = 20MeV and go = 1.38 x lO-^MeV-^. 

FIG. 5. The temperature dependence of the pairing gap. (a) The usual solid metal. We set 
II — m = 5eV and go = 4.42 x 10~^eV~^. (b) The metallic hydrogen. We set ^ — m = 200eV ( 
rs = 0.5 ) and go = 6.90 x 10~^°eV~^ (c) The relativistic plasma. We set /J, — m = 0.5MeV and 
go = 19.7MeV-2. (d) The neutron star. We set fi - m = 20MeV and go = 7.89 x IQ-^MeV-^. All 
calculations are done by taking the cutoff A/{fi — m) = 1. 
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FIG. 6. The presentation of our results for various electron systems in an unified view 
point, by functions of the chemical potential, (a) The coupling constant so as to provide 
|A(r = 0)|/(/Lt — m) = 10^'^ for the solution of the "nonrelativistic" case. By using the data 
given in (a), we calculate (b) the magnitudes of gaps at T = 0, and (c) Tc in 4-cases. All of the 
calculations are performed by taking the cutoff A/{^ — m) = 1. 

FIG. 7. The temperature dependences of spin triplet gap functions, (a) The usual solid metal. 
We set 11 — 171 = 5eV and ^fi = 1.1 x 10~^eV~^. (b) The relativistic plasma. We set — m = 0.5eV 
and gi = 12MeV"^. (c) The neutron star. We set - m = 80MeV and gi = 1.5 x IQ-^MeV"^. 
All calculations are performed with the cutoff A/(/i — m) = 1. The states not shown in this figure 
give extremely small solutions or cannot be solved because of numerical difficulties. 
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TABLES 

TABLE I. Typical order of physical quantities in various systems. The masses are 0.5 MeV for 

the elctron, and 1 GeV for the nucleon. 



system 



particle density (cm ^) ep — m 



|A(r = o)| 



|A(-Z 



usual solid metal^ 
metallic hydrogen*^ 

"Is electron of uranium" gas'^ 
relativistic plasma*^ 
neutron star'' 

^Refs. 7, 8 

''Refs. 9, 10 

^Rcf.l3 

<iRef.ll 

^Ref.l2 



-10^2 

1023 _ 1025 

-10^9 
~ 1032 

1037 _ 10^8 



< lOOK 

~ io2k 



10-3 

10-3 



~ 5eV ~ 10-3 - lO-^eV 
~ 20 - 200eV ~ 10-i - 10-2eV 
~ 0.25MeV 
~ 2.5MeV 

~ 10 - lOOMeV ~ 0.1 - 2.5MeV ~ 0.1 - l.SMeV ~ 10-^ 
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